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Abstract 



Novel method of calculating Nuclear Statistical Equilibrium is presented. 
Basic equations are carefully solved using arbitrary precision arithmetic. Special 
interpolation procedure is then used to retrieve all abundances using tabulated 
results for neutrons and protons, together with basic nuclear data. Proton 
and neutron abundance tables, basic nuclear data and partition functions for 
nuclides used in calculations are provided. Simple interpolation algorithm using 
pre-calculated p and n abundances tabulated as a functions of kT, p and Y e is 
outlined. Unique properties of this method are: (1) ability to pick-up out of 
NSE selected nuclei only (2) computational time scaling linearly with number 
of re-calculated abundances (3) relatively small amount of stored data: only 
two large tables (4) slightly faster than solving NSE equations using traditional 
Newton- Raphson methods for small networks (few tens of species); superior for 
huge (800-3000) networks (5) do not require initial guess; works well on random 
input (6) can tailored to specific application (7) ability to use third-party NSE 
solvers to obtain fully compatible tables (8) encapsulation of the NSE code for 
bug-free calculations. 

Range of applications for this approach is possible: coverage test of tra- 
ditional NSE Newton-Raphson codes, generating starting values, code-to-code 
verification and possible replacement of the old legacy procedures in supernova 
simulations. 
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1. Brief introduction 

Main goal of the article is to provide new method of computing the Nuclear 
Statistical Equilibrium (NSE) abundances of the nuclear species. Vast range 
of conditions can be analyzed: Y e = 0.0 ... 1.0, p = 10 2 . . . 10 13 g/cc and T — 
2 x 10 9 . . . 10 11 K covering almost any astrophysical situation imaginable. While 
this article do not concentrate on particular target object, these results are useful 
for study of pre-supernova stars after Si burning [TT] , thermonuclear supernovae 
[T2] . core-collapse [3], and protoneutron stars pQ. 

We calculate NSE abundances using reliable arbitrary precision arithmetic 
approach. Tables of pre-calculated proton (X p ) and neutron (X n ) abundances 
as functions of the thermodynamic conditions defined by the triad: p, T, Y e are 
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stored. Recovering of the remaining several hundred abundances from these two 
tables is non-trivial task. Detailed description of the working procedure used 
to calculate NSE abundances is provided. Algorithm is fast thanks to use of 
pre-tabulated X pn . It has unique ability to pick up out of NSE ensemble only 
species of interest and other features. Computational time scales linearly with 
number of required nuclides. 



2. NSE 

2.1. Basic equations 

Well-known equations for the ensemble of Ni SO + 1 nuclei in thermal equi- 
librium [HE] are: 

N,„„ 



^2 X k = 1 ( la ) 



fc=0 

JV. 



Y^^X k =Y e (lb) 

fe=0 k 

where abundance X k for fc-th nuclei with atomic number Z k and mass number 
A k is: ^ ^ 

X k = \G k {T) QpAUA 3 ) " A k b l 2 Xt- Zk X^e^. (2) 
Temperature dependent partition function for k-th. nuclei is given by: 

G k (T)=^T(2J lk + l)e-^ (3) 

i=0 

where summation is over all known excited states (numbered by the index i) of 
the k-th nucleus; Ji k and E^ k are the spin and the excitation energy, respectively. 
Q k is the binding energy, p,T - density and temperature of the plasma, Na is the 
Avogadro number and k - Boltzmann constant. Thermal de'Broglie wavelength 
used in eq. ^ is: 

\p2/Kmji kT ^ ^ 

where mn is the mass of the hydrogen atom and h denotes Planck's constant. 

Partition function has been calculated directly from nuclear database using 
pj. Missing spins were assumed to be equal zero. For uncertain lowest possible 
was used. Results are in good agreement with data used in (4j [5] , cf. Fig. [2] 
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Figure 1: Nuclides included in NSE calculations. Much larger (800 nuclear species) networks 
were also tested. 
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Figure 2: Comparison of the partition functions derived directly using ^ and database [B] 
(solid lines), and those from Hix & Thielemann code [1][5] (dashed lines). 



2.2. Limitations of the Newton-Raphson NSE solvers 

NSE equations, from mathematical point of view, form the system of two 
large high order polynomial equations (polynomial system) for unknown pro- 
ton X p and neutron X n abundances. System is solved numerically using two- 
dimensional Newton-Raphson techniqu^] Due to large integer powers and other 
factors this approach is prone to numerous convergence problems. While in "nor- 
mal" situation (typical thermodynamic parameters, good initial guess, standard 
selection of species) convergence of codes using machine floating point arithmetic 
is amazingly fast, failures are inevitable. Limited numerical precision might be 
problematic issue. This forces programmers to include multi-level fail-safe pro- 
cedures. They are by many orders of magnitude slower, and not guaranteed 
to converge. Careful programming with proper handling of round-off errors 
is required to get correct results, leading to additional complications. Due to 
problems with numerical precision and unpredictable iteration numbers, rapidly 
growing with number of species and for low temperatures, procedures are long, 
complicated and hard to parallelize. 

Moreover, even if we are interested in abundance of single nucleus entire 



*In principle polynomial system might be reduced using Groebner basis methods, especially 
over rational field. In practice ensemble including protons, neutrons, 4 He and single heavy 
nuclei can be solved, but additional components cause Groebner basis algorithms to fail in 
sense of computational time: no result is returned in a period of several hours. 
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system ([T]) has to be solved. Such a situation is typical for neutrino spectrum 
calculations, as usually much more nuclear species are included in NSE than 
those with known neutrino emission rates. Usually very few of them contribute 
at non-negligible level, e.g. p, 56 Ni and 55 Co for v e emission at Y e — 0.5. Large 
part of kJ '— p—Y e space is completely dominated by processes involving neutrons 
and protons only. In the course of the research we have faced this problem. In 
recent article [12] NSE ensemble included 800 nuclides while FFN tables used 
include only 189 of them. Interpolation of the pre-calculated results has been 
found to be optimal solution. Similar problem arise in core-collapse supernova 
simulations. Depending on temperature, NSE or full reaction network is solved. 
Again, NSE can be computed for larger ensemble, but due to limitations of 
network ODE solvers only fraction of species is traced. 



2.3. Interpolation algorithm 

To handle results of the NSE calculations efficiently, interpolation seems to 
be wrong solution. Naively, one might try to interpolate stored proton X p and 
neutron X n abundances obtained from Eqns. ([lj, and get X/. from ([2]). Un- 
fortunately, this does not work. Even a very small inaccuracy in X n or X p 
produces enormous error^Jin Xf. due to large (~ A) integer powers in ([2]). An- 
other ,, brute force" method is tabulation of every Xk- This might be useful 
if a few out of NSE species are of interest. This is also the fastest approach. 
However, for larger number of species amount of stored data becomes very large: 
tens or hundreds of tables like Table [7] instead of two. Fortunately, we found 
a compromise, which successfully combines both ideas. Inability to get accu- 
rate abundances using interpolated X n , X p does not include grid points, as they 
can be stored with accuracy equal to the machine precision, or even better if 
required. First, we calculate abundance of selected species Xk at grid points 
neighboring given (p, T, Y e ) point. Next, we interpolate using computed X^s. 
Only proton X p and neutron X p abundances need to be tabulated, but more 
(using formula ^ at 8 corners of a cube) computational time is required com- 
pared to interpolation of stored Xk values for all nuclei. Additionally, partition 
function Gfc(T), atomic and mass numbers Zk,Ak and binding energy Qk has 
to be stored for all nuclei to use pi). Using (tri)linear interpolation eqns. ([!]) 
are fulfilled automatically up to original solving accuracy. 

We still have to solve to generate X p and X p tables. Any method e.g. 
existing codes [TB] , pre-calculated results or a web service [5] may be used in this 
purpose. Because efficiency and speed of the code is not of primary importance 
if one use interpolating scheme, Eqns. (la lb ) has been solved numerically using 
MATHEMATICA code0 Integrated MATHEMATICA \6\ database has been 



2 This relative error can be estimated as: 6(X„ Xp) ~ A2 A SX, where SX is typical 
relative error of X n (X p ) and A is mass number. For A ~ 60 amplification of relative error 
might be as large as 10 18 (!) for X n ~ X p ~ 0.5. 

3 Entire code .9 has approx. 100 lines including database loading, writing C headers, and 
solving Q with arbitrary precision. Code is slow compared to FORTRAN equivalents, a 
price paid for arbitrary precision. This is not an important issue, as all we want is to generate 
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Table 1: Minimum number of nuclides required to compute all abundances above X m i n . 

Z A niso Last included nuclide 
56 Ni 

59 Cu 

60 Zn 

61 Zn 

63 Zn 

63 Ga 

65 Ga 

66 Ge 

68 Ge 
75 Kr 

87 Nb 



28 
28 
29 
30 
30 
30 
31 
31 
32 
32 
36 
41 



56 
57 
59 
GO 
Gl 
G3 
G3 
G5 
GG 
G8 
75 
87 



562 
563 
592 
620 
621 
623 
651 
653 
683 
685 
807 
970 



used, including excited states and spins. This let us to calculate temperature 
dependent partition function. Measured excited states were used if present in 
database, otherwise neglected. Third party partition functions can be used as 
well. No Coulomb and screening corrections were applied. 

Proton and neutron abundances are then tabulated as a functions of temper- 
ature, density and electron fraction. NSE results are checked against available 
codes/results [HI HI [13] with good agreement. 

2-4- Brief discussion of NSE results 

Determination of NSE abundances is crucial for many applications, includ- 
ing nucleosynthesis, neutrino emission, nuclear energy generation and equation 
of state. Therefore we have made some tests to verify results and accuracy es- 
timates. Despite known physical issues (temperature-dependent partition func- 
tion, Coulomb corrections [13], screening [7]) one of the most important factors 
is number and selection of species included in equations (JlJ. Even a single one 
important nuclei missing in NSE ensemble may lead to radically different re- 
sults. While inclusion of some species seems obvious (p, n , 4 He, 56 Ni, iron 
group) further selection is more or less arbitrary. 

To quantify problem I tried to answer the following question: what is the 
maximum required atomic (Z) and mass number (A) to get solution including 
all species with abundance larger than prescribed X m in- Results are presented 
in Table [T] and Fig. [3] For example, from Table [T] if we do not want to miss 
any of species with abundance above e.g. 10 -6 , we need at least nuclides up to 
61 Zn. Nuclei in Fig. [3] are ordered according to [5J; approximate Z and A are 
included as a tick marks for a top axis. This estimate gives an upper limit for 
number of required nuclei. To get true minimal number of nuclides required to 



tables. We do it once, in parallel if required. Later we use interpolators, which are very fast, 
even compared to codes using hardware floats. 
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Figure 3: NSE abundances as a function of the number of the nuclei involved in calculations 
for kT = 0.4 MeV, p = 10 7 g/cm 3 and Y e = 0.5. 



get all species above assumed accuracy one have to consider all subsets for entire 
considered kT — p — Y e space. Number of subsets, given by the Bell number 
Bpf iso is very large. Therefore, rigorous selection of species is impossible for 
large sets, and the safest thing to do is to use estimates given by Table [T] or 
consider all nuclei available [5]. In practice however, other factors decide, e.g. 
limited computational resources in supernova simulations. 

From Fig. [3] we can conclude that the most primitive NSE including p and n 
only is not useful, maybe except for very high temperatures, cf. Fig. [4] Inclusion 
of the alpha particle extends applications to lower temperatures but usually p 
and n abundances are wrong by few orders of magnitude. To get correct abun- 
dances of p and n for lower temperatures entire iron peak has to be included. 
X p and X n are rock-stable if all nuclei below Z=28, A=56 are included. This 
number might be seriously reduced if we focus on narrow Y e range and exclude 
low mass (A=3..16) elements. Anyway, results in Table [l] indicate, that no 
more than 1000 nuclei are required to get all abundances above 10~ 30 . While 
it is possible to solve NSE equation for more than 3000 nuclides [S], it does not 
change results significantly. 
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Ye=0.45, lgp=6 [g/cm 3 ] 




kT [MeV] 

Figure 4: NSE abundance versus temperature. 

2.5. NSE viewgraphs 

We discuss some properties of the NSE state for completeness. For a very 
high temperature^] above kT ~ 0.5 MeV in Fig. [4] (Tg ~5.8) no bound nuclei 
exist and we have a mixture of free neutrons and protons (Fig. [4]). In this case 
solution of system ([I]) is: 

X n — 1 Y e . Xp Y e . 

If temperature decrease helium is being "synthesized" like in Big Bang nucle- 
osynthesis. If temperature drops further more below kT ~ 0.35 MeV (T ~ 
4 x 10 9 K) and thermodynamic conditions are maintained for long enough, 
heavy most bound nuclei are preferred. Finally, cold catalyzed matter state is 
a pure (for Y e = 0.45 ) 58 Fe, cf. Fig. g for Y e = 0.5 it is 56 Ni. This is appeal- 
ing physical picture. Note extremely strong Y e dependence of the NSE state 
(Fig. [5| for 0.35 <Y e < 0.55. Y e dependence for large temperatures is trivial: 
smooth balance between p, n and a abundances. The most interesting is the 
temperature range where heavy nuclei dominate. Note that, for higher densities, 
temperature threshold for heavy nuclei formation moves to higher temperatures, 
see eq. ^ and footnote [i] 



4 Actually, if we forget temperature dependent partition function, according to eq. |5] 
solution depends on factor proportional to p 2 /kT 3 . 
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kT=0.62 MeV, lg p=6 [g/cm 3 ] 




0.0 0.2 0.4 0.6 0.8 1.0 

Figure 5: NSE abundance versus electron fraction Y e . 



Striking feature of Fig. [5] is a rapid variation of the abundances within 
range of Y e = 0.35 . . . 0.55, cf. Fig. [6j NSE prefer nucleus with individual 
Yj k ^ = Z^/Ak as close as possible to Y e for entire thermodynamic ensemble. 
For example double magic nuclei 78 Ni with largest known neutron exces ^(low- 
est Y e — 28/78 ~ 0.36) dominates for Y e < 0.365 (not included in example 
network) until neutrons (with Y e = 0) take a lead. For opposite side, Y e 3> 0.5, 
protons are dominant]^] 

Note, that for Y e = exact solution of NSE equations is X n = 1; for 
Y e = 1 we get X p = 1. Abundances for intermediate values of Y e continuously 
approach these values for Y e — > 0, and Y e — > 1, cf. Fig. [5] Other abundances 
rapidly approach zero. 

Rapid abundance variation has strong imprint on neutrino emission. For ex- 
ample, known for large electron capture rate 55 Co has non- negligible abundance 
only in narrow range of Y e — 0.47 . . . 0.5. 



5 Neutron excess is equivalent to Y e : r\ = 1 — 2Y e . 

6 Normally, for Y e 2> 0.5 protons dominate. But if 3 Li would exist, it should take the 
role of hydrogen under NSE conditions if density is high enough. This species is still present 
in nuclide databases with atomic mass 3.030775 and binding energy 2.2676 MeV, despite 
experimental detection [171 has never been confirmed |15| ; see also comments in ESNDF data 
at http://ie.lbl.gov/ensdf/. 
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kT=0.62 MeV, lgp=6 [g/cm 3 ] 




Figure 6: Zoom of the Fig. [5] into the most interesting range of Y e = 0.35 . . . 0.55. 

3. Proton and neutron NSE abundance tables: explanation and ex- 
amples of use 

Dataset described and presented in the article is meant to be simple exam- 
ple of the methods used. It is tailored to test against 32 isotope NSE solver 
used by Garching group, based on serial code of Hix & Thielemann [H|S]. In 
real application user should use larger tables for bigger networks available on- 
line or generate (request from author) user-defined dataset tailored to specific 
application. 

Here we provide tables of the proton, X p and neutron, X n abundances, 
together with nuclear data required to calculate all remaining abundances, X^. 
Additionally, list of nuclides is required, including: 

1. atomic and mass numbers 

2. masses and binding energies 

3. spins and excited states or, equivalently, tabulated temperature-dependent 
partition function 

In Table [7] we have provided: in the first column temperature kT in MeV, 
second column include base 10 logarithm of the density in g/cm 3 , and third 
column electron fraction Y e , i.e. number of electrons divided by number of 
baryons. Fourth column include proton abundance under NSE defined by the 
kT, p and Y e , and fifth column contain neutron abundance. 
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Table 2: Symbols used in nuclide table, main Table [5] begins on pagc|16| 



Col 


Symbol 


Description 


1 


No 


position 


2 


Symbol 


Standard element symbol 


3 


A 


Mass number 


4 


N 


Neutron number 


5 


Z 


Atomic number 


6 


Q 


Binding energy per nucleon [MeV] 


7 


Jo 


ground state spin (0 if not known) 



Table 3: Symbols used in temperature dependent partition function table, main Table[5]be sins 
on page [17| NOTE: without ground state partition function 2 Jo + 1, included in Table [5] 
Total partition function (J3J is a sum of 2 Jo + 1 and function tabulated below. Results are 
truncated below 10 — 6 . 



Col 


Symbol 


Description 


1 


No 


position 




2 


Symbol 


Standard element symbol 


3 


kT=0.2 


Partition function for kT = 


0.2 MeV 


4 


kT=0.4 


Partition function for kT = 


0.4 MeV 


5 


kT=0.6 


Partition function for kT — 


0.6 MeV 


6 


kT=0.8 


Partition function for kT — 


0.8 MeV 


7 


kT=1.0 


Partition function for kT — 


1.0 MeV 



Table 4: Symbols used in NSE proton and neutron abundance table, main Table [7] begins on 



Col 


Symbol 


Description 


1 


kT 


Temperature [MeV] 


2 


IgP 


base 10 logarithm of the density [g/cm 3 ] 


3 


Y e 


number of electrons per baryon 


4 


X p 


abundance of free protons 


5 


x n 


abundance of free neutrons 
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To calculate all NSE abundance we need basic nuclear data, presented in 
Table [5j and partition functions, from Table [5] 

Using approach presented here, main computational cost is the partition 
function, so use of tabulated version instead of eq. Q is important. 

Detailed description of the algorithm is presented below. Goal is to calculate 
abundance Xk of species k for given temperature T, density p and Y e : 

X k =NSE(T,p,Y e ,k). 

1. from tables of the proton(neutron) abundance we pick up points surround- 
ing requested T,p,Y e ; in case of e.g. trilinear interpolation these points 
are 8 corners of a cuboid - requested point must be inside or at the edge 
of the cuboid 

2. for all these points we calculate abundance Xk from ^ 

3. now we have machine-precision accurate abundances Xi at 8 corners of 
the cuboid 

4. interpolate (trilinear interpolation in the example) to get at desired 
point 

We point out again, that we interpolate X^ NOT X p or X n . X^ must be 
calculated exactly at grid points. Example implementation of the algorithm is 
included in libnse library [TU] , 

3.1. Implementation notes 

Article deals with interpolation of the functions of three variables. Despite 
progress in computer hardware, available memory amount in particular, it is 
hard to find sophisticated 3D interpolators. Therefore trilinear, or mixed bilin- 
ear on T — p plane and staircase for Y e , interpolations were used. 

A lot of computational time is spend on X^X^~ Z . Large integer power of 
the floating-point can be computed nearly optimal using double-exponentiation 
algorithm, usually included in standard math libraries. Minor improvements for 
range of interest can be achieved using optimal chain of powers, C++ template 
programming or other technique devoid of if instruction. 

We still recommend caution with integer powers of floating point numbers. 
For example standard math.h from C does not include integer powers, GNU Gsl 
only up to 9, while cmath.h from C++ standard library does. This cause large 
variation of the computational time. 

Higher-order interpolation might possibly help to fit procedure into CPU 
cache memory thanks to reduced amount of data. However we are also in 
danger of overfitting resulting in catastrophic errors, e.g. negative abundances. 
If amount of memory is not an issue, linear interpolation is recommended. 

Partition function is evaluated using linear interpolation. 

4. Additional numerical data 

Printed tables and results described in the article are meant to be sim- 
ple examples of the proposed method. Extended versions of the tables, custom 
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datasets and numerical library can be downloaded from http://ribes.if.uj.edu. 
or requested from the author. 
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5. NSE proton and neutron abundance table 



Table 5: Nuclei included in NSE and required nuclear data. 



No 


Symbol 


Z 


N 


A 


Q 


Jo 


1 


l n 





1 


1 





1/2 


2 


l U 


1 





1 





1/2 


3 


4 Hc 


2 


2 


4 


7.0739150 





4 


12 C 


6 


6 


12 


7.6801440 





5 


16Q 


8 


8 


16 


7.9762060 





6 


20 Ne 


10 


10 


20 


8.0322400 





7 


24 Mg 


12 


12 


24 


8.2607090 





8 


28 Si 


14 


14 


28 


8.4477440 





9 


32 S 


16 


16 


32 


8.4931340 





10 


36 Ar 


18 


18 


36 


8.5199090 





11 


40 Ca 


20 


20 


40 


8.5513010 





12 


44 Ti 


22 


22 


44 


8.5335180 





13 


50 Ti 


22 


28 


50 


8.7556180 





14 


48 Cr 


24 


24 


48 


8.5722100 





15 


54 Cr 


24 


30 


54 


8.7779140 





16 


55 Cr 


24 


31 


55 


8.7318840 


3/2 


17 


54 Mn 


25 


29 


54 


8.7379230 


3 


18 


55 Mn 


25 


30 


55 


8.7649880 


5/2 


19 


56 Mn 


25 


31 


56 


8.7383000 


3 


20 


52 Fe 


26 


26 


52 


8.6095980 





21 


54 Fe 


26 


28 


54 


8.7363440 





22 


55 Fe 


26 


29 


55 


8.7465600 


3/2 


23 


56 Fe 


26 


30 


56 


8.7903230 





24 


57 Fe 


26 


31 


57 


8.7702490 


1/2 


25 


58 Fe 


26 


32 


58 


8.7922210 





26 


55 Co 


27 


28 


55 


8.6695750 


7/2 


27 


56 Co 


27 


29 


56 


8.6948170 


4 


28 


56 Ni 


28 


28 


56 


8.6427090 





29 


57 Ni 


28 


29 


57 


8.6709010 


3/2 


30 


58 Ni 


28 


30 


58 


8.7320410 





31 


eo Ni 


28 


32 


60 


8.7807570 





32 


60 Zn 


30 


30 


60 


8.5832730 
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Table 6: Temperature dependent partition function 



No 


Name 


0.20 


0.40 


0.60 


0.80 


1.00 


3 


4 He 


0.00 


0.00 


0.00 


0.00 


0.00 


4 


12 C 


0.00 


0.00 


0.00 


0.02 


0.06 


5 


16Q 


0.00 


0.00 


0.00 


0.00 


0.03 


6 


20 Ne 


0.00 


0.08 


0.34 


0.72 


1.21 


7 


24 Mg 


0.00 


0.16 


0.53 


1.00 


1.58 


8 


28 Si 


0.00 


0.06 


0.26 


0.58 


1.00 


9 


32 S 


0.00 


0.02 


0.14 


0.44 


1.02 


10 


36 Ar 


0.00 


0.04 


0.21 


0.61 


1.37 


11 


40 Ca 


0.00 


0.00 


0.04 


0.24 


0.87 


12 


44 Ti 


0.02 


0.38 


1.32 


3.09 


5.90 


13 


50 Ti 


0.00 


0.12 


0.62 


1.79 


4.04 


14 


48 Cr 


0.12 


0.85 


1.98 


3.60 


5.95 


15 


54 Cr 


0.08 


0.74 


2.09 


4.58 


8.62 


16 


55 Cr 


1.39 


5.48 


11.30 


18.60 


27.00 


17 


54 Mn 


10.90 


22.30 


36.70 


54.70 


75.20 


18 


55 Mn 


4.36 


7.81 


14.00 


24.00 


37.40 


19 


56 Mn 


15.80 


33.00 


50.80 


70.20 


90.80 


20 


52 Fe 


0.07 


0.63 


1.49 


2.70 


4.46 


21 


54 Fe 


0.00 


0.19 


1.11 


3.43 


7.57 


22 


55 Fe 


0.33 


2.12 


6.25 


13.20 


22.50 


23 


56p e 


0.07 


0.68 


2.01 


4.78 


9.60 


24 


57 Fe 


7.67 


12.80 


20.20 


30.50 


43.60 


25 


58 Fe 


0.09 


0.87 


2.81 


6.71 


13.10 


26 


55 Co 


0.00 


0.07 


0.74 


2.86 


7.02 


27 


56 Co 


4.07 


10.60 


18.50 


27.90 


38.50 


28 


56 Ni 


0.00 


0.00 


0.08 


0.32 


0.86 


29 


57 Ni 


0.14 


1.04 


2.46 


4.63 


7.87 


30 


58 Ni 


0.00 


0.17 


0.95 


2.93 


6.68 


31 


60 Ni 


0.00 


0.25 


1.31 


3.89 


8.48 


32 


60 Zn 


0.03 


0.46 


1.36 


2.84 


5.03 
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Table 7: Proton and neutron abundance tables. 



kT 




Y e 


X p 


x n 


0.20 


6 





350 


1 


4733451757175688e- 


32 


1.9792286428615166e- 


01 


0.20 


8 





350 


4 


6594595569570535e- 


37 


1.9791728691182955e- 


01 


0.20 


10 





350 


1 


4734610247732251e- 


41 


1.9791672869601915e- 


01 


0.40 


6 





350 


7 


6850895802764179e- 


10 


2.0039781283597258e- 


01 


0.40 


8 





350 


2 


5035658289594687e- 


14 


1.9822605323795905e- 


01 


0.40 


10 





350 


7 


9454130714268810e- 


19 


1.9794846194430299e- 


01 


0.60 


6 





350 


8 


6190577876701473e- 


03 


3.0861905778717863c- 


01 


0.60 


8 





350 


1 


6076815125584695c- 


06 


2.0253972809389420e- 


01 


0.60 


10 





350 


5 


3242575469851928c- 


11 


1.9823081609598006e- 


01 


0.80 


6 





350 


3 


4503551381543635c- 


01 


6.4503551381543633e- 


01 


0.80 


8 





350 


6 


0527987083834160c- 


03 


3.0605116217065886e- 


01 


0.80 


10 





350 


5 


6832537960862639e- 


07 


1.9976073216155885e- 


01 


1.00 


6 





350 


3 


4999839054511411c- 


01 


6.499983905451 1410e- 


01 


1.00 


8 





350 


1 


6548297105954407e- 


01 


4.6548297105954412c- 


01 


1.00 


10 





350 


1 


5943665263060269e- 


04 


2.1299547811356118e- 


01 


0.20 


6 





400 


4 


5278905162835453e- 


32 


8.3343373485753311c- 


02 


0.20 


8 





400 


1 


4321164378624826e- 


36 


8.3334338486210496e- 


02 


0.20 


10 





400 


4 


5288358535831710e- 


41 


8.3333433860011616e- 


02 


0.40 


6 





400 


2 


2541535853182530c- 


09 


8.7066391596373041e- 


02 


0.40 


8 





400 


7 


6454369365453428e- 


14 


8.3829450685805496e- 


02 


0.40 


10 





400 


2 


4404642046330322e- 


18 


8.3384805016069527e- 


02 


0.60 


6 





400 


1 


3274703748466302c- 


02 


2.1327470374824983e- 


01 


0.60 


8 





400 


4 


5119825481078112e- 


06 


9.0877967412558644e- 


02 


0.60 


10 





400 


1 


6254313935731902c- 


10 


8.3849904241914824e- 


02 


0.80 


6 





400 


3 


9448770859198856e- 


01 


5.9448770859198852e- 


01 


0.80 


8 





400 


9 


4260257833551960c- 


03 


2.0942531492521452e- 


01 


0.80 


10 





400 


1 


6763122530124213e- 


06 


8.6551386142692061e- 


02 


1.00 


6 





400 


3 


9999820882445386c- 


01 


5.9999820882445387e- 


01 


1.00 


8 





400 


2 


0024502450742615e- 


01 


4.0024502450742611e- 


01 


1.00 


10 





400 


3 


8713189248855591e- 


04 


1.0652457407487780e- 


01 


0.20 


6 





450 


1 


2254155814326638e- 


15 


1.0461956026905567e- 


14 


0.20 


8 





450 


1 


3863583647470302e- 


17 


1.0941992041239904e- 


16 


0.20 


10 





450 


1 


5862580710493338e- 


19 


1.1335873063097257e- 


18 


0.40 


6 





450 


1 


3854242348476404e- 


05 


8.6981272932959097e- 


05 


0.40 


8 





450 


1 


6020742194861760e- 


07 


8.9771845366122210e- 


07 


0.40 


10 





450 


1 


6200879000788178c- 


09 


1.0323858209007692e- 


08 


0.60 


6 





450 


2 


3969925356851572e- 


02 


1.2396992535682635e- 


01 


0.60 


8 





450 


2 


6113737999714160e- 


04 


3.6480691095702282e- 


03 


0.60 


10 





450 


5 


0716727991727212e- 


06 


2.4807585164625365e- 


05 


0.80 


6 





450 


4 


4414588253856829e- 


01 


5.4414588253856822e- 


01 


0.80 


8 





450 


1 


7653064602620443e- 


02 


1.1765299048431674e- 


01 


0.80 


10 





450 


2 


7450088757524307e- 


04 


1.5448501064190729e- 


03 


1.00 


6 





450 


4 


4999809512767586e- 


01 


5.4999809512767583e- 


01 
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Table 7: Proton and neutron abundance tables [continued]. 



kT 




Y e 


x p 




1.00 


8 


n 
u 




2.4074257595157450e-01 


o 
o 


/in7/19^7^Q^1 £.7/1/1 Sr. 


ni 


1.00 


10 


n 
u 




2.6369076822919937e-03 


o 
z 


oi lUoooooy ( oi / oyc- 


no 


0.20 


6 


n 
u 


ouu 


2.1027565926622363e-06 


q 
o 


3^01 ^0881 fi1 8QGG8/3 


ZO 


0.20 


8 





500 


2.0902933564115093c-07 


3 


9833138940328108e- 


26 


0.20 


10 





500 


2.1502284101490003e-08 


4 


5645516281556141c- 


29 


0.40 


6 





500 


2.4656257836791810e-02 


1 


1948909248170071e- 


07 


0.40 


8 





500 


2.9898277397867462e-03 


1 


2462253458017922e- 


10 


0.40 


10 





500 


2.7224923203959757e-04 


1 


6226416771 168784e- 


13 


0.60 


6 





500 


5.5105593212499610e-02 


5 


5105593212499276e- 


02 


0.60 


8 





500 


2.7482483740179343e-02 


6 


7168188134800034e- 


05 


0.60 


10 





500 


6.4480794989162076e-03 


3 


9682763842094195e- 


08 


0.80 


6 





500 


4.9402975852718162e-01 


4 


9402975852718162e- 


01 


0.80 


8 





500 


4.6130285406649901e-02 


4 


6130284779868531e- 


02 


0.80 


10 





500 


2.1479785529572089e-02 


3 


4028673055878073e- 


05 


1.00 


6 





500 


4.9999805645135403e-01 


4 


9999805645135403e- 


01 


1.00 


8 





500 


2.8751054543310173e-01 


2 


8751054543310173c- 


01 


1.00 


10 





500 


3.0258682045056728e-02 


2 


9362486076929912e- 


03 
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